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Abstract 

The criticality of the (2 + l)-dimensional XY model is investigated with the numerical diagonal- 
ization method. So far, it has been considered that the diagonalization method would not be very 
suitable for analyzing the criticality in large dimensions (d > 3); in fact, the tractable system size 
with the diagonalization method is severely restricted. In this paper, we employ Novotny's method, 
which enables us to treat a variety of system sizes N = 6, 8, . . . , 20 (TV: the number of spins con- 
stituting a cluster). For that purpose, we develop an off-diagonal version of Novotny's method to 
adopt the off-diagonal (quantum-mechanical XY) interaction. Moreover, in order to improve the 
finite-size-scaling behavior, we tune the coupling-constant parameters to a scale-invariant point. 
As a result, we estimate the critical indices as v = 0.675(20) and 7/1/ = 1.97(10). 
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I. INTRODUCTION 



It has been considered that the diagonalization method would not be very suitable for 
analyzing the criticality in large dimensions d > 3. In fact, as the system size enlarges, the 
number of spins constituting a cluster increases rapidly in d > 3, and the dimensionality of 
the Hilbert space soon exceeds the limitation of available computer resources. Such a severe 
limitation as to the tractable system size prevents us from making a systematic analysis of 
the simulation data. 

To cope with this difficulty, Novotny proposed a transfer- matrix formalism [l, 2, 3], which 
enables us to construct a transfer-matrix unit with an arbitrary (integral) number of spins 
N; note that conventionally, the number of spins is restricted within N = 2 d ~ 1 , 3 d_1 , .... As 
a demonstration, Novotny simulated the Ising model in d < 7 systematically 3j. Meanwhile, 
it has been shown that the idea is applicable to a wide class of systems such as the frustrated 
!sing model |4] and the quantum-mechanical Ising model under the transverse magnetic field 

A 

In this paper, we extend the Novotny method to adopt the off-diagonal (quantum- 
mechanical XY) interaction; see the Hamiltonian, Eq. ([T]), mentioned afterward. Actually, 
as mentioned above, the use of Novotny's method has been restricted within the case of the 
diagonal (Ising-type) interaction. As a demonstration, we apply the method to the (2 + 1)- 
dimensional XY model with a variety of system sizes N = 6, 8, . . . , 20. Taking the advantage 
that a series of system sizes are available, we made a systematic finite-size-scaling analysis 
of the simulation data. As a result, we estimate the critical indices as v = 0.675(20) and 
7/V = 1.97(10). Recent developments on the d = 3 XY universality class are overviewed 
with an emphasis on the microgravity-environment experiment; see also Refs. 




11]. Our method provides an alternative approach to the d = 3 XY universality 



To be specific, we consider the following Hamiltonian for the (2 + l)-dimensional XY 



model 



12 



13 



14j with the extended interactions 



H 



J \ fqxqx,qyqy\ J 
-<JNN 3j) — J NNN 



(ij) {(ij)) [ijkl] i 

(1) 

Here, the quantum spin-1 (S = 1) operators {Si} are placed at each square-lattice point 



z\2 
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i. The summations, an d ^2\ijkih run over au possible nearest-neighbor, next- 

nearest-neighbor, and plaquette spins, respectively. The parameters, Jnn, Jnnn, and D n , 
are the corresponding coupling constants. The single- ion anisotropy D, drives the system 
from the XY phase (D < D c ) to the large- .D phase (D > D c ). (In the large- .D phase, the 
ground state is magnetically disordered, accompanied with a finite excitation gap.) Our aim 
is to survey the criticality by means of the off-diagonal Novotny method. 

The Hamiltonian ([1]) has a number of tunable parameters. We fixed them to 

(Jnn, Jnnn, D n ) = (0.158242810160, 0.058561393564, 0.10035104389), (2) 

and survey the D-driven phase transition. As explicated in the Appendix, around the point 
([2]), the finite-size-scaling behavior improves significantly; the irrelevant interactions cancel 
out, because the point (T2J) is a scale- invariant point with respect to the real-space decimation 
shown in Fig. [91 Such an elimination of finite-size corrections has been utilized successfully 



to analyze the criticality of the classical systems such as the Ising model [15|, [16] and the 
lattice 4 theory 17, [l8] . We adopt the idea to investigate a quantum-mechanical system, 
Eq. ©. 

The rest of this paper is organized as follows. In Sec. [Til we develop an off-diagonal 
version of the Novotny method. In Sec. IIII[ employing this method, we simulate the (2 + 1)- 
dimensional XY model (CQ). In Sec. IIV[ we present the summary and discussions. In the 
Appendix, we determine a scale-invariant point (T5]) with respect to the real-space decimation 
shown in Fig. [9l 

II. OFF-DIAGONAL NOVOTNY'S METHOD 

In this section, we explain the simulation scheme. As mentioned in the Introduction, we 
develop an off-diagonal version of Novotny's method to simulate the (2 + l)-dimensional XY 
model (ED) ; so far, the Novotny method has been applied to the case of Ising-type interactions 

To begin with, we explain the basic idea of Novotny's method. Novotny's method allows 
us to construct a cluster with an arbitrary number of spins; see Fig. [TJ As indicated, the 
basic structure of the cluster is one- dimensional. The dimensionality is lifted to d = 2 by the 
bridges over the (V / A)-th-neighbor pairs. Because the basic structure is one-dimensional, 
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we are able to construct a cluster with an arbitrary (integral) number of spins {Si} (i = 
1,2, ... , N); note that naively, the number of spins is restricted within N — 4, 9, 16, . . . in 
d = 2. 

We formulate the above idea explicitly We propose the following expression 

N 

n = -j nn (h xy (i)+h xy (Vn))-j nnn (h xy (Vn+i)+h(Vn-i))+^^ 

8=1 

(3) 

for the Hamiltonian of the (2 + l)-dimensional XY model ([T]). The component Hxy(u)(v) 
describes the XY (plaquette) interaction over the v-th neighbor pairs; see Fig. [TJ Because 
the quantum XY interaction, Hxy{v), is an off-diagonal one, we need to develop an off- 
diagonal version of Novotny's method. We propose the following expression 

N 

H X y{v) = J2( pVS i P ~ VS i + P v SfP~ v Sf). (4) 

i=l 

This formula serves as a basis of the off-diagonal Novotny method. The symbol P denotes 
the translation operator by one lattice spacing; 

P\Si, 5*2, . . . , Sn) = \Sn, Si, ... , Sjv-i)- (5) 

(We impose the periodic boundary condition, SV+i = Si.) Here, the base |{Sfc}) diagonalizes 
the {S-} operators; namely, it satisfies 

S?\{S h }) = St\{S k }), (6) 

for each I = 1, 2, . . . , N. The insertions of the operators P in Eq. (j3J) introduce the 
v-th neighbor interaction along the alignment of spins {Si}; symbolically, the operator 
P v SfP~ v may be written as S^ +v . On the other hand, as for H n (v), we adopt the con- 
ventional idea based on the diagonal Novotny method That is, its diagonal elements 
{({S k }\H n (v)\{S k })} are given by 

({Sk}\Ha(v)\{S k }) = ({S k }\P°T\{S k }), (7) 
with the four-spin interaction 

N 

({S k }\T\{T k }) = J2 S ^+^ T ^- ( 8 ) 
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Similarly, the insertion of P v introduces the v-th neighbor interaction. However, in the 
diagonal scheme ([7j), one operation of P v suffices; note that in the off-diagonal formalism (j4j), 
two operations P are required. Because each operation P ±v requires huge computational 
effort, the off-diagonal scheme is computationally demanding. Afterward, we provide a 
number of formulae useful for the practical implementation of the algorithm. 

The above formulae complete the formal basis of our simulation scheme. Aiming to 
improve the simulation result, we implement the following symmetrization technique 
That is, we symmetrize the component Hxy,u{v) by replacing it with 

H X Y,n(v) -> (H XY:D (v) + H XY ,u(-v))/2. (9) 

This replacement restores the symmetry between the ascending, Si, 5*2, . . . , SV, and the 
descending, Sat, Sn-i, ■ ■ ■ , Si, directions completely. 

Last, we provide a number of formulae that may be useful in the practical implementation 
of the algorithm. We utilize the translationally invariant bases {\k,n)}, which diagonalize 
the operator P; 

P\k,n) = e ik \k,n). (fO) 

Here, the wave number k runs over a Brillouin zone k = 2nM/N (M: integer), and the 
index n specifies the state within the subspace k. As anticipated, the bases {\k,n)} are 
useful to obtain an explicit representation of the formulae mentioned above. For instance, 
the first term of the formula (j3J) is represented by 

Af Af 

(k,n\J2 S j pVS j P ~ V \ k > m ) = ^^{k ) n\S x ] \k' 1 n'){k' 1 n'\S x j \k ) m)e i{k '- k)v 1 (ff) 

j=l j=l k',n' 

in terms of the frame {\k,n)}. Because the parameter v is, in general, an irrational num- 
ber, the oscillating factor e^ k '~ k ' v is incommensurate with respect to the lattice periodicity. 
Hence, the intermediate summation Y2k' nas t° ^ e treated carefully; namely, each Brillouin 
zone {k'} is no longer equivalent. We accepted the following symmetrized sum 

o-k' — 2 a ~ n a -K+2Tr/N + a -7r+47i7Af + ■ • • + o-tt-2-k/n + 2 a ?>"' (12) 

k' 

with respect to a summand ay. Here, the denominators of the first and the last terms 
compensate the duplicated sum at the edges of the Brillouin zone [— (Similarly, we 
obtain an explicit representation for H n (v) via the conventional Novotny method l|, |2j.) 



Provided that the explicit matrix elements of Hxy,d{v) are at hand, we are able to perform 
the numerical diagonalization of the Hamiltonian ([3]). The results are shown in the next 
section. 

Last, we make an overview of the S — 1 XY model (TXT) . As mentioned in the Introduc- 

nnn 

tion, the model has been studied in Refs. [12j, [13j, [14|]. In the case of d = 1 dimension, the 
criticality (D-driven phase transition) was investigated in detail According to Ref. 19]. 
for sufficiently large D, a magnetically disordered ground state (large- .D phase) appears, and 
the criticality is identical to that of the classical counterpart in d + 1(= 2) dimensions (KT 
transition). Unfortunately, a naive extension to the case of S = 1/2 is not appropriate, 
because the D anisotropy, D(S*) 2 , reduces to a constant term, D/4. (Moreover, the trans- 
verse magnetic field violates the XY symmetry, and the criticality changes into the Ising 
type.) As a matter of fact, it is difficult to realize a ground-state phase transition for the 
S = 1/2 model without violating the translational invariance and the rotational symmetry. 
(Possibly, the double-plane S = 1/2 model may exhibit a desirable criticality by tuning the 
inter-plane interaction. However, this model is too complicated.) Hence, we consider the 
S = 1 XY model (CQ) with the .D-anisotropy term. 



III. NUMERICAL RESULTS 



In this section, we analyze the criticality of the (2 + l)-dimensional XY model, Eq. (CD). 
As mentioned in the Introduction, the coupling-constant parameters ( Jnn, Jnnn, D n ) are 
set to the scale-invariant point ([2]). Thereby, we survey the D-driven phase transition with 
the finite-size scaling. In order to diagonalize the Hamiltonian, we utilize the off-diagonal 
Novotny method developed in Sec. [Til Owing to this method, we treat a variety of system 
sizes AT = 6, 8, ... , 20. The linear dimension L of the cluster is given by 

L = VN, (13) 

because the N spins constitute a rectangular cluster; see Fig. [TJ 



A. Transition point 

In this section, we provide an evidence of the D-driven phase transition, and estimate 
the critical point D c with the finite-size scaling. 
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In Fig. [21 we plot the scaled energy gap LAE for various D, and N = 6, 8, . . . , 20 with 
the other coupling constants fixed to Eq. ([2]). The symbol AE denotes the first-excitation 
gap. According to the finite-size scaling, the scaled energy gap LAE should be invariant 
at the critical point. Indeed, we observe an onset of the .D-driven phase transition around 
D« 1. 

In Fig. [31 we plot the approximate transition point D c {Li ) L 2 ) for (2/(Li + L2)) 3 with 



6 < N\ < N 2 < 20 and Li >2 = -JN\^\ the validity of the l/L 3 -extrapolation scheme (abscissa 
scale) is considered at the end of this section. Here, the approximate transition point 
D c (Li,L 2 ) denotes a scale-invariant point with respect to a pair of system sizes (Li,L 2 ). 
Namely, the approximate transition point satisfies the equation 

LiAE{L\)\d = d c {l u L2) = L 2 AE(L 2 )\ D=Dc(LlM) . (14) 

The least-squares fit to the data of Fig. [3]yields D c = 0.9569(83) in the thermodynamic limit, 
L — > 00. As a reference, we calculated D c = 0.9744(68) through the l/L 4 -extrapolation 
scheme. Considering the deviation as an error indicator, we estimate the critical point as 

D c = 0.957(25). (15) 

Let us mention a few remarks. First, we consider the abscissa scale 1/L 3 utilized in 
Fig. [3j Naively, the scaling theory predicts that dominant corrections to D c should scale 
like with uo = 0.785(20) and v = 0.6717(1) [11]. O n one hand, in our simulation, 

such dominant corrections should be suppressed by tuning the coupling constants to Eq. 
[2]) ; see the Appendix. The convergence to the thermodynamic limit may be accelerated 
16f| . (For extremely large system sizes, the singularity 1/L W+1 ^ may emerge.) Hence, in 
Fig. [31 we set the abscissa scale to 1/L 3 . Second, we argue a consistency between the 
finite-size scaling and the real-space decimation; in the Appendix, we made a fixed-point 
analysis (lAllj) . regarding D as a unit of energy D = 1 (1A7I) . This proposition D = 1 is 



quite consistent with the above scaling result (fl~5T) . validating the fixed-point analysis in the 
Appendix. In other words, around the fixed point (1A11I) . corrections to scaling may cancel 
out satisfactorily. Encouraged by this consistency, in Sec. IHI C[ we survey the criticality 
rather in detail. 
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B. Comparison with the conventional XY model 

In the preceding section, we simulated the XY model ([1]) with the finely tuned coupling 
constants As a comparison, in this section, we provide the data for the conventional 
XY model. That is, we turn off the extended coupling constants, setting the interactions 
to (Jnn, Jnnn, Dn) = (0.2,0,0) tentatively. 

In Fig. HI we plot the scaled energy gap LAE for various D and N — 6,8, ... , 20. We 
observe an onset of the D-driven phase transition around D m 1.1. However, the data are 
scattered, as compared to those of Fig. [2j In fact, in Fig. HI the location of the transition 
point appears to be less clear. This result demonstrates that the finely-tuned coupling 
constants ([2]) lead to elimination of finite-size corrections. 



C. Critical exponents 

In Sec. IHI A[ we observe an onset of the D-driven phase transition. In this section, we 
calculate the critical exponents, v and 7/V, based on the finite-size scaling. 
In Fig. [51 we plot the approximate critical exponent 

v(T t) = HLi/L 2 ) 

1 h 2) \n(d D (L 1 AE(L 1 ))/d D (L 2 AE(L 2 )))\ D=D : { ) 

for (2/(Li + L 2 )f with 6 < N x < N 2 < 20 (L lj2 = y/N^, and D c = 0.957 [Eq. 05]]; 
afterward, we consider the abscissa scale, 1/L 2 . The least-squares fit to these data yields 
v = 0.675(16). As a reference, we calculated v = 0.687(11) through the l/L 3 -extrapolation 
scheme. Considering the deviation as an error indicator, we estimate the critical exponent 

as 

v = 0.675(20). (17) 
In Fig. El we plot the approximate critical exponent 

1 /u = \n(x±(L 1 )/x±(L 2 ))\ D=Dc /HL 1 /L 2 ), (18) 

for (2/(L! + L 2 )) 2 with 6 < N x < N 2 < 20 (L lj2 = y/N^), and D c = 0.957 [Eq. JEJ]. Here, 
the transverse susceptibility x± is given by the resolvent form 

Xx = ^(alMxj^-^MM, (19) 
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with the ground state \g) and the ground-state energy E g . The magnetization M x is given 
by M x = £^ =1 Sf. The resolvent form (fl9j) is readily calculated with use of the continued- 



fraction method 



20|. 



The least-squares fit to the data in Fig. [6] yields j/i/ — 1.965(61). As a reference, 
we calculated 7/z/ = 2.020(42) through the l/L 3 -extrapolation scheme. Considering the 
deviation as an error indicator, we estimate the critical exponent as 

7/z/ = 1.97(10). (20) 

Last, we argue the abscissa scale 1/L 2 utilized in Figs. [5] and [61 Naively, the scaling 
theory predicts that dominant corrections to the critical indices should scale like 1/ L u with 
uj = 0.785(20) [11]. On one hand, as argued in Sec. IIII Al such dominant corrections 
should be suppressed by adjusting the coupling constants to Eq. ([2]), and the convergence 



is accelerated than the naively expected one 16j. Hence, we set the abscissa scale to 1/L 2 
in Figs. [5] and [61 



D. Refined data analysis 

In this section, we make an alternative analysis of the criticality to demonstrate a relia- 
bility of our scheme. 

In Figs. [7J and [HI we plot the critical exponents 



v 



(21) 



D=D C (L U L 2 ) 

and 

7/z/ = Hxl{Li)Ixx.{L2))\d=d c{ l„l 2 )IHLiIL 2 ), (22) 

respectively, for (2/ (Li + L 2 )) 2 with 6 < Ni < N 2 < 20. Here, these exponents are calculated 
at the approximate critical point D = D c (Li,L 2 ) (THl) rather than at D c = 0.957 as in the 
preceding section. 

Clearly, these data, Figs. [7] and [HI exhibit accelerated convergence to the thermodynamic 
limit as compared to those of Figs. [5] and [61 In fact, the least-squares fit to these data yields 
the estimates v = 0.658(5) and 7/z/ = 1.946(4) with suppressed error margins. Actually, 
in Fig. [HI the systematic error dominates the insystematic one. In such a case, one has to 
make a detailed consideration of the nature of corrections to scaling to ensure the accuracy 



(amount of error margin) of the extrapolation. Here, we do not commence making such a 
consideration, and accept the estimates, Eqs. ffTT|) and fl20l . obtained less ambiguously in 
the preceding section. It is not the purpose of this paper to obtain fully refined estimates 
for the critical indices. Such a detailed analysis will be pursued in the succeeding works. In 
fact, the diagonalization method has a potential applicability to the frustrated magnetism, 
for which the quantum Monte Carlo method suffers from the notorious sign problem. The 
Novotny method would be particularly of use to explore such a problem. Actually, in the 



case of the Ising-type anisotropy, the Novotny method was applied [2l| to clarifying the 
nature of the frustration-driven transition (Lifshitz point). The present scheme may provide 
a basis for surveying such a quantum frustrated system with the AY-type symmetry. 

IV. SUMMARY AND DISCUSSIONS 

The criticality of the (2 + 1) -dimensional XY model (0Q) was investigated with the 
numerical-diagonalization method. For that purpose, we developed an off-diagonal ver- 
sion of Novotny's diagonalization method (Sec. [Til), which enables us to treat a variety of 
system sizes N = 6,8, ...,20 (N: the number of spins within a cluster). Moreover, we 
improved the finite-size-scaling behavior by adjusting the coupling-constant parameters to 
a scale-invariant point (j5J). 

Owing to these improvements, we could analyze the simulation data systematically with 
the finite-size scaling. As a result, we estimated the critical indices as v = 0.675(20) and 
7/1/ = 1.97(10). These indices immediately yield the following critical exponents 

a = -0.025(60), (5 = 0.348(49), and 7 = 1.330(78), (23) 

through the scaling relations. 

Recent developments on the d = 3 XY universality class are overviewed in Ref. 
Our diagonalization result ( |23l) is accordant with a Monte Carlo result, a = —0.0151(3) 
f3 = 0.3486(1), and 7 = 1.3178(2) [ll|, and a field-theoretical result, a = -0.011(4), (3 = 



0.3470(16), and 7 = 1.3169(20) [8]. (In Ref. Ill ], information from a series-expansion result 
is also taken into account.) To the best of our knowledge, no numerical-diagona ization 
result has been reported as for the d = 3 XY universality class. According to Ref. jg], there 
arose a discrepancy between the Monte Carlo simulation and the microgravity-environment 
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experiment; see also Ref. [23j. As a matter of fact, the microgravity experiment [22J reports a 
critical exponent a = —0.0127(3). In order to resolve this discrepancy, an alternative scheme 
other than the Monte Carlo and series-expansion methods would be desirable. Refinement 
of the diagonalization scheme through considering the singularity of corrections to scaling 
might be significant in order to settle this longstanding issue. 
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APPENDIX A: SEARCH FOR A SCALE-INVARIANT POINT: ELIMINATION 
OF FINITE-SIZE CORRECTIONS 

As mentioned in the Introduction, we simulated the quantum XY model ([IT), setting 
the coupling constants to Eq. (jSJ); around this point, we observe eliminated finite-size 
corrections. In this Appendix, we explicate the scheme to determine the point (TSJ). 

To be gin with, we explain the technique to suppress the finite-size corrections. According 

m 



to Rcfs. 11!: 



17 



18] , the finite-size behavior improves around the renormalization-group 
fixed point. That is, the irrelevant interactions may cancel out around the fixed point. 
Clearly, such an improvement of the finite-size behavior admits us to make a systematic 
finite-size-scaling analysis of the simulation data. To avoid confusion, we stress that the 
fixed-point analysis is simply a preliminary one, and subsequently, we perform large-scale 
computer simulation to estimate the critical exponents. In this sense, as for the Monte 
Carlo simulation, it might be more rewarding to enlarge the system size rather than to 
extend the coupling-constant parameters and adjust them. On one hand, it is significant 
for the numerical diagonalization to eliminate corrections to scaling, because its tractable 
system size is restricted intrinsically. 

In Fig. M, we present a schematic drawing of the real-space-decimation procedure. As 
indicated, we consider a couple of rectangular clusters with the sizes 2x2 and 4x4. These 
clusters are labeled by the symbols 5* and L, respectively. Decimating out the spin variables 
indicated by the symbol • within the L cluster, we obtain a coarse-grained lattice identical to 
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the S cluster. Our aim is to search for a scale-invariant point with respect to the real-space 
decimation. 

Before going into the fixed-point analysis, we set up the simulation scheme for the clusters, 
S and L. We cast the Hamiltonian ([T]) into the following plaquette-based expression 

-h = J2^m + dJ2(sD 2 , (Ai) 

[ijkl] i 

with the plaquette interaction 

JNN / q x q x qy qy q x q x . qy qy , qx QX i QV QV I QX QX i QV Q1J\ ( \ 0\ 

n ijkl — 2 ^ « j °i °j °i °/ °j °/ "r 2 fe °i °fcJ l^" Z J 

-Jnnn(S?S? + SfSf + Sp x k + Spl) (A3) 
+J D D (^ + ^ + ^ + ^) 2 . (A4) 

(The denominator of the coefficient Jnn compensates the duplicated sum.) Hence, the 
Hamiltonian for the S cluster is given by 

4 

-Hs = H° 2U + dY j (S:) 2 , (A5) 

i=l 

with the replacement 

J N N ,N N N — > (1 + b)JNN,NNN- (A6) 

Here, the parameter 6 controls the boundary interaction strength, and hereafter, we set 
b = 0.7; we consider the validity of this choice afterward. The boundary-interaction param- 
eter b interpolates smoothly the open, 6 = 0, and periodic, 6 = 1, boundary conditions. The 
point is that for the two-site (L = 2) system, the bulk interaction, SfS%, and the boundary 
interaction S^S" coincide each other. Hence, for the S cluster, the boundary interaction b is 
freely tunable without violating the translation invariance. We make use of this redundancy 
to obtain the fixed point reliably. On the other hand, the L cluster does not have such a 
redundancy, and the Hamiltonian TIl is given by Eq. ([T|) with L = 4 unambiguously. We 
diagonalize these Hamiltonian matrices TCs,l numerically; note that we employ the conven- 
tional diagonalization method, rather than the off-diagonal Novotny method developed in 
Sec. H 

With use of the simulation technique developed above, we search for the fixed point of 
the real-space decimation. We survey the parameter space (Jnn, Jnnn, D n ), regarding D 
as a unit of energy; namely, we set 

D = 1, (A7) 
12 



throughout this section. Thereby, we impose the following conditions 



2AE S = AAE L (A8) 
(5*5f) s = (§Z§$)l (A9) 
(5-51)5 = (A10) 

as a scale-invariance criterion. The symbol AE$,l denotes the excitation gap for the respec- 
tive clusters. The arrangement of the spin variables, Sf 2 3 4 an d Sf 2 3 a, is shown in Fig. [91 
The symbol (. . .)s,l denotes the ground-state average for the respective clusters. The first 
equality (IA8I) comes from the scale invariance of the scaled energy gap, LAE. (We refer 



readers to Ref. 



241 ]. where the author utilizes such a critical-amplitude relation success- 



fully to analyze the renormalization-group flow numerically.) On one hand, the remaining 



equations, (1A9j) and (1A10j) . are the scale-invariance conditions 25j regarding the correlation 



functions for the edge and diagonal spins, respectively. 

The conditions, Eqs. (1A8I) - (1A10I) . are the nonlinear equations with respect to 
( Jnn, Jnnn, D n ). In order to obtain the solution, we employed the Newton method, and 
found that the following nontrivial solution does exist; 

(Jnn, Jnnn, D a ) = (0.158242810160, 0.058561393564, 0.10035104389). (All) 

The last digits may be uncertain because of the round-off errors. 

Last, we argue the validity of the above solution (lAllj) and the boundary condition 



b = 0.7. In Sec. 1111} via the finite-size-scaling analysis, we obtained D c = 0.957(25) (fT5l) . 
Apparently, this result is consistent with D = 1 postulated in Eq. (1A7I) . Moreover, the 
simulation data in Fig. [2] exhibit suppressed finite-size corrections, as compared to those 
of the ordinary XY model, Fig. [61 These features validate the choice of the boundary 
condition b = 0.7 as well as the reliability of the fixed point (lAllj) . Furthermore, we point 



out that the boundary condition b = 0.7 is reminiscent of b = 0.4 utilized in the fixed-point 
analysis of the d = 3 Ising ferromagnet . 
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FIG. 1: A schematic drawing of the spin cluster for the d = 2 quantum XY model ([I]) is presented. 
As indicated above, the spins constitute a one-dimensional alignment {Si} (i = 1,2, ... , N), and 
the dimensionality is lifted to d = 2 by the bridges over the (7V 1 /2) t h -neighbor pairs. This is a 
basic idea of Novotny's diagonalization method. We need to develop an off-diagonal version of 
Novotny's method, because we have to adopt the quantum XY interaction ([!]); see Sec. HH 




FIG. 2: Scaled energy gap LAE is plotted for various D and N = 6, 8, . . . , 20 (L = \N); note that 
we survey the D-driven phase transition with the other interactions, (Jnn, Jnnn, D n ), adjusted 
to a fixed point ([2]). We observe a clear indication of the D-driven transition around D ~ 1. 
Apparently, the finite-size-scaling behavior is improved as compared to that of the conventional 
XY model (Fig. 
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FIG. 3: The approximate critical point D c {UD is plotted for (2/(Li+L 2 )) 3 with 6 < Ni < N 2 < 20 
(Li t 2 = y/Nip). The least-squares fit to these data yields D c = 0.9569(83) in the thermodynamic 
limit L — > cxd. 
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FIG. 4: Tentatively, we turned off the extended interactions, (Jnn , Jnnn , D a ) = (0.2,0,0), and 
calculated the scaled energy gap LAE for various D and N = 6, 8, . . . , 20 (L = VN). We notice 
that the data are scattered as compared to those of Fig. [2l 
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FIG. 5: The approximate critical exponent v (|16p is plotted for (2/ '(L1+L2)) 2 with 6 < N% < N2 < 
20 (Li,2 = \JN\j2). The least-squares fit to these data yields v = 0.675(16) in the thermodynamic 
limit L — > 00. 
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FIG. 6: The approximate critical exponent 7/V (I18h is plotted for (2/{L\ + L2)) 2 with 6 < 
N\ < N2 < 20 (-^1,2 = y/Nip)- The least-squares fit to these data yields "f/v = 1.965(61) in the 
thermodynamic limit L — » 00. 
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FIG. 7: The approximate critical exponent v (|21|) is plotted for (2/ '(L1+L2)) 2 with 6 < N\ < N 2 < 
20 (^1,2 = \J~Nip)- The least-squares fit to these data yields v = 0.658(5) in the thermodynamic 
limit L — > 00. 
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FIG. 8: The approximate critical exponent 7/V (1221) is plotted for {2/{L\ + L2)) 2 with 6 < 
N\ < N 2 < 20 (Li )2 = y/N lt2 ). The least-squares fit to these data yields j/v = 1.946(4) in the 
thermodynamic limit L — » 00. 
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FIG. 9: A schematic drawing of the real-space renormalization group (decimation) for the d = 2 
XY model ([1]) is presented. Through decimating out the spin variables indicated by the symbol 
• within the L cluster, we obtain a coarse-grained lattice identical to the S cluster. Imposing the 
scale-invariance conditions, Eqs. (|A8p - (|A10p . we arrive at the fixed point ([2]). 
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